161 research outputs found
Corner and finger formation in Hele--Shaw flow with kinetic undercooling regularisation
We examine the effect of a kinetic undercooling condition on the evolution of
a free boundary in Hele--Shaw flow, in both bubble and channel geometries. We
present analytical and numerical evidence that the bubble boundary is unstable
and may develop one or more corners in finite time, for both expansion and
contraction cases. This loss of regularity is interesting because it occurs
regardless of whether the less viscous fluid is displacing the more viscous
fluid, or vice versa. We show that small contracting bubbles are described to
leading order by a well-studied geometric flow rule. Exact solutions to this
asymptotic problem continue past the corner formation until the bubble
contracts to a point as a slit in the limit. Lastly, we consider the evolving
boundary with kinetic undercooling in a Saffman--Taylor channel geometry. The
boundary may either form corners in finite time, or evolve to a single long
finger travelling at constant speed, depending on the strength of kinetic
undercooling. We demonstrate these two different behaviours numerically. For
the travelling finger, we present results of a numerical solution method
similar to that used to demonstrate the selection of discrete fingers by
surface tension. With kinetic undercooling, a continuum of corner-free
travelling fingers exists for any finger width above a critical value, which
goes to zero as the kinetic undercooling vanishes. We have not been able to
compute the discrete family of analytic solutions, predicted by previous
asymptotic analysis, because the numerical scheme cannot distinguish between
solutions characterised by analytic fingers and those which are corner-free but
non-analytic
Efficient computation of two-dimensional steady free-surface flows
We consider a family of steady free-surface flow problems in two dimensions,
concentrating on the effect of nonlinearity on the train of gravity waves that
appear downstream of a disturbance. By exploiting standard complex variable
techniques, these problems are formulated in terms of a coupled system of
Bernoulli's equation and an integral equation. When applying a numerical
collocation scheme, the Jacobian for the system is dense, as the integral
equation forces each of the algebraic equations to depend on each of the
unknowns. We present here a strategy for overcoming this challenge, which leads
to a numerical scheme that is much more efficient than what is normally
employed for these types of problems, allowing for many more grid points over
the free surface. In particular, we provide a simple recipe for constructing a
sparse approximation to the Jacobian that is used as a preconditioner in a
Jacobian-free Newton-Krylov method for solving the nonlinear system. We use
this approach to compute numerical results for a variety of prototype problems
including flows past pressure distributions, a surface-piercing object and
bottom topographies.Comment: 20 pages, 13 figures, under revie
Spectrograms of ship wakes: identifying linear and nonlinear wave signals
A spectrogram is a useful way of using short-time discrete Fourier transforms
to visualise surface height measurements taken of ship wakes in real world
conditions. For a steadily moving ship that leaves behind small-amplitude
waves, the spectrogram is known to have two clear linear components, a
sliding-frequency mode caused by the divergent waves and a constant-frequency
mode for the transverse waves. However, recent observations of high speed ferry
data have identified additional components of the spectrograms that are not yet
explained. We use computer simulations of linear and nonlinear ship wave
patterns and apply time-frequency analysis to generate spectrograms for an
idealised ship. We clarify the role of the linear dispersion relation and ship
speed on the two linear components. We use a simple weakly nonlinear theory to
identify higher order effects in a spectrogram and, while the high speed ferry
data is very noisy, we propose that certain additional features in the
experimental data are caused by nonlinearity. Finally, we provide a possible
explanation for a further discrepancy between the high speed ferry spectrograms
and linear theory by accounting for ship acceleration.Comment: 21 pages, 10 figures, submitte
Numerical investigation of controlling interfacial instabilities in non-standard Hele-Shaw configurations
Viscous fingering experiments in Hele-Shaw cells lead to striking pattern
formations which have been the subject of intense focus among the physics and
applied mathematics community for many years. In recent times, much attention
has been devoted to devising strategies for controlling such patterns and
reducing the growth of the interfacial fingers. We continue this research by
reporting on numerical simulations, based on the level set method, of a
generalised Hele-Shaw model for which the geometry of the Hele-Shaw cell is
altered. First, we investigate how imposing constant and time-dependent
injection rates in a Hele-Shaw cell that is either standard, tapered or
rotating can be used to reduce the development of viscous fingering when an
inviscid fluid is injected into a viscous fluid over a finite time period. We
perform a series of numerical experiments comparing the effectiveness of each
strategy to determine how these non-standard Hele-Shaw configurations influence
the morphological features of the inviscid-viscous fluid interface. Tapering
plates in either converging or diverging directions leads to reduced metrics of
viscous fingering at the final time when compared to the standard parallel
configuration, especially with carefully chosen injection rates; for the
rotating plate case, the effect is even more dramatic, with sufficiently large
rotation rates completely stabilising the interface. Next, we illustrate how
the number of non-splitting fingers can be controlled by injecting the inviscid
fluid at a time-dependent rate while increasing the gap between the plates.
Simulations compare well with previous experimental results for various
injection rates and geometric configurations. Further, we demonstrate how the
fully nonlinear dynamics of the problem affect the number of fingers that
emerge and how well this number agrees with predictions from linear stability
analysis
A sharp-front moving boundary model for malignant invasion
We analyse a novel mathematical model of malignant invasion which takes the
form of a two-phase moving boundary problem describing the invasion of a
population of malignant cells into a population of background tissue, such as
skin. Cells in both populations undergo diffusive migration and logistic
proliferation. The interface between the two populations moves according to a
two-phase Stefan condition. Unlike many reaction-diffusion models of malignant
invasion, the moving boundary model explicitly describes the motion of the
sharp front between the cancer and surrounding tissues without needing to
introduce degenerate nonlinear diffusion. Numerical simulations suggest the
model gives rise to very interesting travelling wave solutions that move with
speed , and the model supports both malignant invasion and malignant
retreat, where the travelling wave can move in either the positive or negative
-directions. Unlike the well-studied Fisher-Kolmogorov and Porous-Fisher
models where travelling waves move with a minimum wave speed ,
the moving boundary model leads to travelling wave solutions with . We interpret these travelling wave solutions in the phase plane and
show that they are associated with several features of the classical
Fisher-Kolmogorov phase plane that are often disregarded as being nonphysical.
We show, numerically, that the phase plane analysis compares well with long
time solutions from the full partial differential equation model as well as
providing accurate perturbation approximations for the shape of the travelling
waves.Comment: 48 pages, 21 figure
Saffman-Taylor fingers with kinetic undercooling
The mathematical model of a steadily propagating Saffman-Taylor finger in a
Hele-Shaw channel has applications to two-dimensional interacting streamer
discharges which are aligned in a periodic array. In the streamer context, the
relevant regularisation on the interface is not provided by surface tension,
but instead has been postulated to involve a mechanism equivalent to kinetic
undercooling, which acts to penalise high velocities and prevent blow-up of the
unregularised solution. Previous asymptotic results for the Hele-Shaw finger
problem with kinetic undercooling suggest that for a given value of the kinetic
undercooling parameter, there is a discrete set of possible finger shapes, each
analytic at the nose and occupying a different fraction of the channel width.
In the limit in which the kinetic undercooling parameter vanishes, the fraction
for each family approaches 1/2, suggesting that this 'selection' of 1/2 by
kinetic undercooling is qualitatively similar to the well-known analogue with
surface tension. We treat the numerical problem of computing these
Saffman-Taylor fingers with kinetic undercooling, which turns out to be more
subtle than the analogue with surface tension, since kinetic undercooling
permits finger shapes which are corner-free but not analytic. We provide
numerical evidence for the selection mechanism by setting up a problem with
both kinetic undercooling and surface tension, and numerically taking the limit
that the surface tension vanishes.Comment: 10 pages, 6 figures, accepted for publication by Physical Review
Time-frequency analysis of ship wave patterns in shallow water: modelling and experiments
A spectrogram of a ship wake is a heat map that visualises the time-dependent
frequency spectrum of surface height measurements taken at a single point as
the ship travels by. Spectrograms are easy to compute and, if properly
interpreted, have the potential to provide crucial information about various
properties of the ship in question. Here we use geometrical arguments and
analysis of an idealised mathematical model to identify features of
spectrograms, concentrating on the effects of a finite-depth channel. Our
results depend heavily on whether the flow regime is subcritical or
supercritical. To support our theoretical predictions, we compare with data
taken from experiments we conducted in a model test basin using a variety of
realistic ship hulls. Finally, we note that vessels with a high aspect ratio
appear to produce spectrogram data that contains periodic patterns. We can
reproduce this behaviour in our mathematical model by using a so-called
two-point wavemaker. These results highlight the role of wave interference
effects in spectrograms of ship wakes.Comment: 14 pages, 7 figure
The role of initial geometry in experimental models of wound closing
Wound healing assays are commonly used to study how populations of cells,
initialised on a two-dimensional surface, act to close an artificial wound
space. While real wounds have different shapes, standard wound healing assays
often deal with just one simple wound shape, and it is unclear whether varying
the wound shape might impact how we interpret results from these experiments.
In this work, we describe a new kind of wound healing assay, called a sticker
assay, that allows us to examine the role of wound shape in a series of wound
healing assays performed with fibroblast cells. In particular, we show how to
use the sticker assay to examine wound healing with square, circular and
triangular shaped wounds. We take a standard approach and report measurements
of the size of the wound as a function of time. This shows that the rate of
wound closure depends on the initial wound shape. This result is interesting
because the only aspect of the assay that we change is the initial wound shape,
and the reason for the different rate of wound closure is unclear. To provide
more insight into the experimental observations we describe our results
quantitatively by calibrating a mathematical model, describing the relevant
transport phenomena, to match our experimental data. Overall, our results
suggest that the rates of cell motility and cell proliferation from different
initial wound shapes are approximately the same, implying that the differences
we observe in the wound closure rate are consistent with a fairly typical
mathematical model of wound healing. Our results imply that parameter estimates
obtained from an experiment performed with one particular wound shape could be
used to describe an experiment performed with a different shape. This
fundamental result is important because this assumption is often invoked, but
never tested
- …